Nonadiabatic coherent evolution of two-level systems under spontaneous decay 
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Abstract 

In this paper we extend current perspectives in engineering reservoirs by producing a time-dependent master equation leading 
to a nonstationary superposition equilibrium state that can be nonadiabatically controlled by the system-reservoir parameters. 
Working with an ion trapped inside a nonindeal cavity we first engineer effective Hamiltonians that couple the electronic states 
of the ion with the cavity mode. Subsequently, two classes of decoherence-free evolution of the superposition of the ground and 
decaying excited levels are achieved: those with time-dependent azimuthal or polar angle. As an application, we generalise the 
purpose of an earlier study [Phys. Rev. Lett. 96, 150403 (2006)], showing how to observe the geometric phases acquired by 
the protected nonstationary states even under a nonadiabatic evolution. 



In the last decade, research activity on the subject 
of open quantum systems has been mainly devoted to 
the search for mechanisms to bypass decoherence. Be- 
yond the quest for conditions that weaken the system- 
reservoir coupling pj, coherence control schemes have 
been introduced among the protocols for quantum error- 
correcting codes @, the existence of decoherence-free 
subspaces (DFS) in collective systems and dynam- 
ical decoupling (DD) methods (4j. More recently, a tech- 
nique has been presented to reach the same goal as the 
DD schemes, without interfering directly in the system 
within the reservoir timescale [5[. We finally mention 
the engineering reservoir program [6], where a quantum 
system whose state is to be protected is compelled to 
engage in additional interactions besides that with the 
reservoir. Such interactions are carefully engineered to 
modify the Liouvillian in a specific way that drives the 
system to an equilibrium with the reservoir. The engi- 
neering reservoir has been developed for both trapped 
ions 0, Q and atomic two- level (TL) systems We 
stress that the engineering reservoir is deeply connected 
to the engineering Hamiltonian program which has been 
pursued for quantum information purposes [lfj. More 
recently, under the assumption of a squeezed engineered 
reservoir, a way to observe the adiabatic geometric phase 
acquired by a protected state evolving coherently through 
the adiabatic manipulation of the squeeze parameters of 
the engineered reservoir has been proposed [ll|, [l2|, [l3| . 

In this paper we show how to protect a nonstation- 
ary superposition state, broadening the range of the pro- 
posed scheme for engineering reservoirs [6J. We present 
a general recipe to build nonadiabatic coherent evolu- 
tions driven by engineered reservoirs. Given this general 
recipe, the task to achieve a particular nonadiabatic evo- 



lution — through a particular engineered reservoir - 
relics entirely on the engineering Hamiltonian program. 
Differently from the standard application of the engineer- 
ing reservoir technique [H, 0, Q , in our model we achieve 
the nonadiabatic decoherence-free evolution of superposi- 
tion states, and show how to implement it in a particular 
system. Working with an ion trapped inside a nonideal 
cavity, we initially engineer an effective Hamiltonian cou- 
pling the electronic levels of the ion with the cavity mode. 
Assuming a strong decay rate of the cavity field, this ef- 
fective interaction is employed to build an artificial reser- 
voir, leading to two classes of asymptotic nonstationary 
superpositions of the ground \g) and decaying excited |e) 
ionic levels: those with time-dependent azimuthal and 
with time-dependent polar angles. We stress that by 
combining both evolutions we can perform any trajec- 
tory on the Bloch sphere. As an application of this en- 
gineering reservoir technique, we generalize the protocol 
in Ref . [n], [3| , demonstrating how to observe geometric 
phases acquired by protected nonstationary states even 
under a nonadiabatic evolution. As a particular case of 
our model, we also show how to build a quantum mem- 
ory to protect stationary superpositions of the internal 
degrees of freedom of the ion. 

The main goal of the standard engineering reservoir 
scheme [6] is to obtain, in the interaction picture, a mas- 
ter equation in the form 



p = - (20p0t 



(1) 



where T is the effective decay rate of the engineered reser- 
voir which is coupled to the quantum system in a spe- 
cific way characterized by the time-independent system 
operator O. The only pure steady state of this system 
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is the eigenstate \ip) of the operator O with null eigen- 
value, ensuring that there is no further eigenstate \<j>) 
of O such that [0,0^] \<j)) =0 0. Even considering 
time-dependent Lindblad operators in master equation 
([1]), the scheme used to protect a given state remains 
valid, assuming an adiabatic evolution of the reservoir 
parameters, i. e., the rate of change of operator O, char- 
acterized by w, is much smaller than T. Consequently, 
we obtain a nonstationary protected state \ip (t)), which 
is the instantaneous eigenstate of O with null eigen- 
value and follows the adiabatic changes of the reservoir 
parameters [l4| • Of course, the fidelity of the protected 
state in this adiabatic evolution depends on the ratio 
zu/T <C 1. Next, to remove the adiabatic constraint in 
the decoherence-free evolution described above, we con- 
sider the engineered time-dependent master equation in 
the interaction picture (H = 1) 

p = -i[H (t) ,p} + ^ [20 (t) pO^ (i) - 0t {t) o (t) p 

- P OHt)0(t)], (2) 

where the Hermitian Hamiltonian H (t) must be 
chosen in accordance with the time dependence of 
the operator O (t) = R (t) OrR) (t), with R(t) = 

Texp^— i f£ H (t 1 ) dt'^j , T being the time-ordering op- 
erator. Note that through the unitary transformation 
R (t) , we recover the time- independent form of the mas- 
ter equation given in (fTJ), in a representation where O 
is replaced by Or. Interestingly, the protected station- 
ary eigenstate \ipn) (Or \iPr) — 0), turns out to be a 
nonstationary state in the original interaction picture, 
\ip(t)) = R(t) \iPr), whose evolution can be manipulated 
by means of appropriate engineered Hamiltonian H(t) 
and reservoir. 

Now we show how to implement the ideas discussed 
above, using a TL trapped ion characterized by the tran- 
sition frequency loq between the ground \g) and excited 
|e) states and trap frequency v. The transition \g) <-» |e) 
is driven by (one or two) classical fields of frequencies u>e, 
with coupling strengths fig (with £ = 1,2), and the ion is 
made to interact — under the Jaynes-Cummings Hamil- 
tonian and Rabi frequency g — with a cavity mode of 
frequency tu a . Within the rotating- wave approximation 
(RWA), the Hamiltonian modelling the system is given 
by 



a a) a + LUQa z /2 + vb^b + |g cos ( k .~x^j 



H = uj 



^ gi(fei- x +01— Ult) ^i(fc 2 . 2; +02— "2*J 



H.c. 



(3) 

where a) (a) and b^ (b) are the creation (annihilation) 
operators of the quantized harmonic field and the vi- 
brational mode whose position operator is x = (b^ + 
b)/\/2mvx 1 m being the ionic mass. The wave vec- 
tors k , k 1 , and k 2 stand for the cavity mode and 



the two amplification fields (with dephasings 4>i and 
</> 2 ), respectively, while au = \k) (l\ (k and I being the 
states g and e). The vibrational mode is decoupled from 
the remaining degrees of freedom of our model by as- 
suming the wave vectors k , k 1 , and k 2 to be per- 
pendicular to He. Under this assumption and going 
to the interaction picture through the transformation 
U = exp [— i (u! a a^a + ujQU z j2^ t\ , we end up with the 
transformed Hamiltonian 



;ff +H.c. 
(4) 

where we have defined the detunings A a = uj a — luq and 

A e — LUi - LOq. 

Nonadiabatic decoherence-free evolution . To accom- 
plish a decoherence-free evolution of a superposition of 
the atomic levels, we must first engineer the appro- 
priate interaction between these levels and the cavity 
mode. To this end we have to adjust the first classi- 
cal field to resonance with the atomic transition, i.e., 
Ai = 0. In what follows, we perform two consecutive 
unitary transformations, first to a framework rotating 
with frequency fli, U\ = exp [— (e 1 ^ 1 cr eg + H.c.) t] , 
which is straightforwardly done with the help of the 
basis states (|±) = (|e) ± e - ^ 1 \g)) /\/2}, constituting 
the eigenstates of the operator Six (e*^ 1 a eg + H.c.) 
defining U±. The adjustment A2 = — 2f2i en- 
ables us to proceed to the second transforma- 
tion U 2 = exp [-ifi 2 (e lv 0-+- + H.c.) t/2] , per- 
formed with the help of another set of basis states 

= (|+) ± e~ lv |— )) /v^j, composed by eigenstates 

of the operator (e 1 ^ +H.c). Working in the 

regime where f2i ~ A 2 ^> f2 2 ~ A a ^> g and adjusting 
A a = —Vt2 : we end up, after a RWA, with the effective 
Hamiltonian 



Ho 



(5/2) 



atcr U 



acr ltj 



(5) 



Now, with the engineered interaction (0 and the dis- 
sipative mechanisms of both the harmonic mode and the 
TL system, the evolution of the transformed density op- 
erator p of the whole system is given by 



p = -i [H 2 ,p] + (r/2) {2apa) - a)ap - pa) a) + CtlP, 

where CtlP stands for the Liouvillian dynamics of the 
TL system under the transformations which do not 
modify the usual Liouvillian form for the harmonic field 
decay. Towards the engineered reservoir, we assume that 
the decay constant of the cavity field is significantly larger 
than both the effective coupling g/2 and the decay con- 
stant 7 of the TL system in CtlP- In our "cavity QED 
+ trapped ion" system, the regime T ^ g, 7 is eas- 
ily achieved through a cavity with low quality factor 
Q = Lu a /F- Together with the good approximation of 
a reservoir at absolute zero, this regime enables us to 
consider only the matrix elements p mn = (m\ p \n) inside 
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the subspace {|0),|1)} of photon numbers. Moreover, 
following the reasoning in Ref . [8j , the strong decay rate 
r enables the adiabatic elimination of the elements poi 
and pn, prompting the evolution of the TL system 

Ptl = r e „ 9 (2(7 Ti p Ti crx T - @ iiPtl ~ Ptl&u) + CtlPtl, 

(6) 

where T eng = g 2 /T stands for the coupling strength of 
the engineered reservoir. The inevitable and undesired 
action of the multimode vacuum CtlP thus works against 
the protected state |f) (f| of the TL system which follows 
asymptotically from Eq. ([6]) with 7 = 0. In fact, taking 
into account the multimode vacuum, the equations of 
motion of the matrix elements (ptl) ts = (f \ P \s) (with r 
and s standing for f and J.), following from the rotating- 
wave approximation, are given by 



(ptl) n = ( T en g + 37/8) - (T eng + 67/T) (ptl) tt 
(ptl) = - ( W2 + 5 7 /4) ( P tl) u + (7/8) (ptl) 



,r 



Ptl 



IT 



whose asymptotic solution leads to the protected state 



p ri (*-foo) = (l-e)|T> <TI+e||> (|| . 



(7) 



where e = [2 + (8/3) (T eng /j)] \ From Eq. ©, it is 
straightforward to compute the fidelity of the protected 
state |T), given by T = Tr[||) (T|prz(* oo)] = 1 - 
e, which approaches unity for a successfully engineered 
coupling strength T eng 3> 7. It is worth noting that even 
the modest ratio T eng /j = 10 results in a fidelity around 
97%. Considering that the coupling between the ground 
and excited states is induced by a Raman transition Q, 
where g ~ 10 5 s" 1 and 7 « 10 2 s _1 HHH, we obtain for a 
cavity decay constant T ss 10 6 s^ 1 the strength T eng /^ « 
10 2 . Therefore, with the excellent approximation e <C 1 
and reversing the unitary transformations U2 and U\ , we 
note that the state |f), written in the interaction picture 
as 

\ip(t)) = cos (ip/2 - flit) \e) - i e-^ 1 sin(^/2 - tlrf) \g) , 

(8) 

allows for a nonadiabatic coherent evolution of a TL sys- 
tem under spontaneous decay, which can be manipulated 
through the laser parameters fii, tp — fa — fa. Such 
an evolution corresponds to flips in the atomic states, 
representing trajectories on different meridian planes on 
the Bloch sphere, governed by the unitary evolution 
exp{-i [(fii*- ip/2) ((Tege^ 1 +H.c.)]}. At this point 
we note that the time-independent operator an corre- 
sponds to Or mentioned above, while the transformation 
U2U1 corresponds to R. 

Quantum memory and adiabatic decoherence-free 
evolution. To build a quantum memory, a device 



that protects a stationary superposition, we turn 
off the second amplification field (f2 2 = 0) in our 
starting Hamiltonian (J4j) - Applying the unitary 
transformation U± = exp [«Ai<j z /2t], we thus obtain 
the time-dependent Hamiltonian Hi = Ai<j z /2 + 
fii (e* o~ eg + a ge ) + [g e^ iA °* aa eg + H.c] . 
Under the additional transformation U 2 = 
exp {-1 [Aitr z /2 + Cli (e* a eg + e"* a ge )] t}, ac- 
complished with the help of the basis states 
{\±) = (y/I^\e)±e-**WJrx\g))/2}, with 
X = Ai/A, A = + Af /4, and the detuning 

A a = — 2A, we finally obtain the effective interaction 



H 2 = (g/2) (e* aV+_ + e~ #1 a<7_ + ) 



(9) 



where lj = g(\ — x/2). 

Following the steps outlined above for the addi- 
tion of the damping mechanism to the cavity mode 
and the TL system, we reach the master equa- 
tion p TL = r e „g (2a + ^p TL a- + - a — Ptl - Ptlo — ) + 
CtlPtl-i with T eng = g^/T, leading to the asymptotic 
solution 

2r £ (i-oo) = (l-e)|+) (+|+r[-><-| 



where e 



e(l-i)- 
1- ff ens /7j 



<C 1 for the case where 



Teng 3> 7, providing again a fidelity ( T = 1 — e) around 
unity for the protected state |+). For the approximation 
(s <C 1), the state |+) written in the interaction picture 
is given by 



In the resonant case, where Ai = x = 0, we obtain the 
stationary state |+) = (|e) + e -1 ^ 1 \g)) /v2, which would 
have been its complementary state |— ) if we had set the 
detuning A a = 2A or the phase fa — > fa + tt. In that 
case, it is worth noting that if the value of the dephasings 
changes adiabatically between fa and fa + 2tt, in such a 
way that the system is always in equilibrium with the 
engineered reservoir, the Bloch vector representing the 
protected state (fTU|) performs a complete rotation around 
the Bloch sphere, as required in Ref. [ll[, to achieve a 
coherent evolution of a superposition state driven by an 
engineered reservoir. 

From the above definition, we note that the ratio 
Teng/l = — x/2)] 2 / (jT) is a function of the param- 
eter x which defines the polar angle of the state vector 
(fTT)|) on the Bloch sphere. This ratio reaches a maximum 
when Qx/A.i — > with negative Ai and, consequently, 
X = — 2, corresponding to the ground state \g). When 
fii/Ai — > with positive Ai, such that x — 2, the 
ratio Teng/7 is null, forbidding us from protecting the 
excited state |e). Evidently, the values Ai,x ^ de- 
scribes nonadiabatic evolution on the Bloch sphere. For 
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the intermediate value x — Oj we obtain from the typi- 
cal strengths considered above for g, 7, and T the value 
Tengh = 10 2 , representing a fidelity around 99% for 
the protected 'equatorial' state |+). Whereas the 'equa- 
torial' case Ai = x — can De employed to achieve 
an adiabatic azimuthal evolution of the stationary state 
1+) = (l e ) + c - ^ 1 \g)) /\/2 on the Bloch sphere, the cases 
Ai,X 7^ describe nonadiabatic evolution on the parallel 
planes on the Bloch sphere. 

Geometric phase induced by reservoir. As can 
be seen from Eqs. ([5]) and (fit))) , in the context 
of decoherence-free evolution, we are able to engi- 
neer the non-stationary superposition of atomic states, 
evolving coherently and acquiring geometric and dy- 
namic phases. In fact, rewriting the state ([8} as 
|^(t)) = [j+> + e -i(v-2Ui*) |_)] /V2, we obtain after 
a cyclic evolution (T = n/fti) the dynamic phase 
4>d(T) = - / T (V(t)| Hj{t) ^ = -tt0 2 / (2fii) and 

the geometric one 0g c '(T) = i Jq (^{t)\ j- t \ip(t)) dt = 
—it, respectively, where Hj{t) = VL\ (0++ — cr__) + 
^ 2 [e 1 ^- 201 *)^- + e - 4 ^- 2f2lt )o-_+] /2. Therefore, 
considering the total evolution time T, under the regime 
of parameters stated above (Q.2 <C we see that the 
contribution coming from the dynamic phase is negligi- 
ble, while the geometric phase is it. 

In order to observe geometric effects we consider an 
auxiliary atomic level a, which does not couple with 
the states \g) and |e) through the action of the fields 
involved in the engineering scheme. We observe that 
within the time scale T of the experiment, the lifetime 
of state I a) does not affect the dynamics described by 
the master equation ©. Otherwise the level a can be 
chosen as a more excited metastable state. To measure 
the phases acquired by the state \ip(t)), we must em- 
ploy an interferometric scheme with \a) as the reference 
state [ll|. For this purpose, using the conservation of 
the total probability, pi a {t) + p+ + {t) + pL-(t) = 1, we 
solve the system of coupled differential equations for the 
probability amplitudes pfj(t) (i,j = +,—,a) following 
from Eq. © in the interaction picture, disregarding, 
within the evolution time T, the small contribution of 
Ctl- Supposing now that the initial state of the sys- 



tem is |*/(0)) = (|^(0)) + |a))/\/2,with 0i = 02 = 0, 
we find at time t that p 1 {t) = | where 
|4r z (f)) = (\a) + e-KHi-KVa)* ^(^j As we have 

set cf>i = (p2 = 0, the protected state at t — turns out 
to be \ip(0)) — \e). For this reason, the superposition 
state 1^/(0)) may be obtained by applying a laser pulse 
between the states \a) and |e). Finally, the geometric 
phase may be observed through the population inversion 
P ea (t) = cos [(2fii + fi 2 ) t] /2 ~ cos [20g c -(t)] /2, where 
t = rm/Clv,n £ N. 

It is worth noting that, differently from the scheme 
proposed in Ref. where the superposition (t)), 

used to measure the geometric phase, is affected by de- 
coherence under a nonadiabatic evolution, here (t)} is 
unaffected by the reservoir even under such an evolution 
faster than that determined by the time scale of the en- 
gineered reservoir. 

Summarizing, we have improved the engineering reser- 
voir by producing a time-dependent master equation 
leading to a nonstationary superposition equilibrium 
state that can be nonadiabatically controlled by the 
system-reservoir parameters. Working with an ion 
trapped inside a bad cavity we constructed two classes 
of decoherence-free evolution of the ground and excited 
ionic levels. By combining the two classes of evolution, 
we can manipulate trajectories on the Bloch sphere by 
changing, alternately, the polar and azimuthal angles. 
Although in our schemes the protected states acquire 
dynamic phases, this fact is unimportant, since they re- 
main in a decoherence-free subspace, where the dynamic 
and geometric phases are unaffected by the reservoir. Fi- 
nally, we have also generalized the objective of the Refs. 
[HE!) showing how to observe the geometric phases ac- 
quired by the protected nonstationary states even under a 
nonadiabatic evolution. We believe that the extention of 
the present scheme for the nonadiabatic time-dependent 
control of a set of qubits, generating quantum logic op- 
erations inside decoherence-free subspaces, may improve 
quantum computation. 
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